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Amorphous solids that underwent a strain in one direction such that they responded in a plastic 
manner 'remember' that direction also when relaxed back to a state with zero mean stress. We 
address the question 'what is the order parameter that is responsible for this memory?' and is 
therefore the reason for the different subsequent responses of the material to strains in different 
directions. We identify such an order parameter which is readily measurable, we discuss its trajectory 
along the stress-strain curve, and propose that it and its probability distribution function must form 
a necessary component of a theory of elasto-plasticity. 



I. INTRODUCTION 

An amorphous solid which is freshly produced by cool- 
ing a glass-forming system from high to low temperature 
is isotropic up to small statistical fluctuations. Put un- 
der an external strain, its stress vs. strain curve should 
exhibit symmetry for positive or negative strains. This 
is not the case for the same amorphous solid after it 
had been already strained such that its stress exceeded 
its yield-stress where plastic deformations become nu- 
merous, resulting in an elasto-plastic flow state. The 
phenomenon is clearly exhibited in Fig. 1. A typical 
averaged stress-strain curve for a 2-dimensional model 
amorphous solid (see below for numerical details) start- 
ing from an ensemble of freshly prepared homogenous 
states is shown in the left panel, with a symmetric tra- 
jectory for positive or negative shear strain. Once in 
the steady flow state, each system in the ensemble is 
brought back to a zero-stress state, which serves as the 
starting point for a second experiment in which a posi- 
tive and negative strain is put on the system as shown 
in the right panel of Fig. 1. Even though the initial en- 
semble is prepared to have zero mean stress, the aver- 
age trajectory now is asymmetric, with positive strain 
exhibiting 'strain hardening' [1], but reaching the same 
level of steady state flow-stress, whereas, the negative 
strain results in a 'strain softening' and a faster yield 
with eventually reaching the same value of steady-state 
flow-stress (in absolute value). This simple phenomenon, 
sometime referred to as the Bauschinger effect [2] , shows 
that the starting point 70 for the second experiment (re- 
ferred below as the Bauschinger point) retains a memory 
of the loading history, some form of anisotropy, which is 
the subject of this article. We stress that the issue under 
study is different from anisotropic elasticity which stems 
from, say, a lattice anisotropy of a crystalline solid. Here 
the systems under study are amorphous, and neverthe- 
less develop a strain induced anisotropy which is much 
more subtle to identify and quantify. 

How to identify the order parameter which is responsi- 
ble for the anisotropy underlying the Bauschinger effect 
is a question that hovers in the elasto-plastic community 
for some while [3-5]. One obvious concept, i.e. of 'back 
stress' [6] or 'remnant stress' for explaining the asymme- 
try seen in the second experiment in Fig. 1 can be ruled 
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FIG. 1: Color online: Stress-strain curves. Left panel: start- 
ing the experiment from a freshly prepared sample results in 
a symmetric trajectory for 7 — ► —7. Right panel: starting the 
experiment from the zero-stress state with 7 = 70 results in 
an asymmetric trajectory, see text for details. Data was av- 
eraged over 500 independent stress-strain curves at T = 0.01 
where temperature is measured in units of e/fes, see Subsect 
II B below. 

out simply by verifying that the initial point has zero 
mean stress. A more sophisticated proposition is em- 
bodied in the 'shear-transformation zone' theory (STZ) 
in which it is conjectured that plasticity occurs in local- 
ized regions whose densities differ for positive and nega- 
tive strains, denoted n + and n_ [7, 8]. The normalized 
difference between these, denoted as m, is a function of 
the loading history and can, in principle, characterize the 
anisotropy that we are seeking. Unfortunately the precise 
nature of the STZ's was never clarified, and it is unknown 
how to measure either n+, n— or m, making it quite im- 
possible to put this proposition under a direct test. More 
recently it was proposed that the sought after anisotropy 
can be characterized in granular matter by the fabric ten- 
sor F = (nn) which captures the mean orientation of the 
contact normals, n, through the spatial average of their 
diadic product [9] . This order parameter was generalized 
for silica glass where n was chosen as a unit vector in the 
direction of the vector distance between Si atoms, disre- 
garding the oxygens. Attempting to test this proposition 
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in the context of the best-studied model of glass-forming, 
i.e. a binary mixtures of point particles with two inter- 
action lengths, or in the case of multi-dispersed point 
particles (see below for details), did not reveal any sys- 
tematic signature of anisotropy. We thus conclude that 
this order parameter is not sufficiently general to be of 
universal use in the development of the theory of elasto- 
plasticity, and that the question of identifying a missing 
order parameter remains open. 
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FIG. 2: Color Online: Model stress-strain curves as obtained 
from cr ~ croc tanh(/i|7 — 7o|/cToo) + P\l - 7o| 2 e~' 7 ~ 70 with 
/3 positive. 



II. THE PROPOSED ORDER PARAMETER 

To see what may serve as a general order parameter we 
examine first the situation with the isotropic amorphous 
solid which is obtained after a quench without any load- 
ing history. Denoting the shear stress by a and the shear 
strain by 7, we observe that isotropy dictates that all 
the even derivatives d 2n a j <hj 2n must vanish by symme- 
try. For example a function that can model the stress- 
strain curve with this constraint in mind may be a ~ 
(Too tanh(/^7/cr 00 ) where a ^ is the flow stress (the mean 
value of the stress in the elasto-plastic steady state) , and 
\i is the shear modulus. For 7 — > —7 this function is per- 
fectly anti-symmetric as required for an isotropic system. 
Imagine now that we add even derivatives to this func- 
tion, say a ~ Coo tanh(/i7/cr 00 ) + /?7 2 e~ 7 with j3 having 
the dimension of stress. The effect will be to change the 
stress-strain curve as seen in Fig. 2, which is quite rem- 
iniscent of the Bauschinger effect. We therefore propose 
that it is advantageous to focus on the even derivatives 
of a vs. 7, with the most important one being the second 
derivative. Note nevertheless that the mechanism leading 
to the existence of a second derivative in our systems is 
not obvious in this simple model. The second derivative 
appears due to plastic deformations whose effect adds 
up to breaking the isotropic symmetry of the freshly 
quenched state. We will show that at the Bauschinger 



point the second derivative is non-zero due to existing 
closer mechanical instabilities in one straining direction 
than in the opposite. 

A. Statistical Mechanics 

Under external loads the displacement field v describes 
how a material point moved from its equilibrium position. 
The strain field is defined (to second order) as 

2 \dxp dx a dx a dxpj 

where here and below repeated Greek indices are summed 
upon. We expand the free energy density J-/V up to a 
constant in terms of the strain tensor 

(2) 

The mean stress is defined as a a p = y ~§r~~^ > an d 

In our simulations we apply a simple shear deformation 
using the transformation of coordinates according to 

Xi -> Xi + 5-iyi , 

Hi — > yi , [%i — > Zi in 3 dimensions) 

where ^7 = 7 — 7* is a small strain increment from any 
reference strain 7*. The explicit 2D strain tensor follow- 
ing Eq. (1) is 

- 1 f 5 "t \ (K\ 

with an obvious generalization in 3D. Since e xx = 0, the 
mean shear stress reduces to the form (equally valid in 
2D and 3D) 

a xy = C* y + C* yxy S 1 +±{C* yyy + Cl yxyxy )5 1 2 + 0(6 1 3 ) . 

(6) 

As discussed above, in isotropic systems where <j xy is 
antisymmetric in Sj, Cf = and the sum C^ yyy 4- 
g X y X y X y _ q Qui proposition is to use the athermal 
limit of this sum as the characterization of the anisotropy 
that we seek. 

B. Models and numerical procedures 

Below we employ a model system with point parti- 
cles of equal mass m and positions in two and three- 
dimensions, interacting via a pairwise potential of the 
form 




where is the distance between particle i and j, e is the 
energy scale, and x c is the dimensionless length for which 
the potential will vanish continuously up to q derivatives. 
The coefficients C2t are given by 
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We chose the parameters k — fO, q — 2 and x c — 1 .385. 
In the three-dimensional simulations each particle i is 
assigned an interaction parameter Ai from a normal dis- 
tribution with mean (A) and Ay = j(Ai + Xj). The 
variance is governed by the poly-dispersity parameter 



A = 15% where A 2 



In the two dimen- 



sional simulations we use the same potential but choose 
a binary mixture model with 'large' and 'small' particles 
such that \ LL = 1.4, X LS = 1.18 and Ass = 1-00. Below 
the units of length are A = Ass m 2D and A = (A) in 
3D. We measure energy, mass and temperature in units 
of e, m and {e/ks) respectively. The unit of time is 
t = yfm\ 2 /e. In the 3D simulations below the mass 
density p = mN/V = 1.3, whereas in 2D p = 0.85. 
In all cases the boundary conditions are periodic and 
thermostating is achieved with the Berendsen scheme 
[12]. We employ the shod equations of motion for im- 
posing deformations, and integrate them using a stan- 
dard leap-frog algorithm [12]. The strain rate is cho- 
sen to be 7 = 10 -4 t _1 for all simulations described be- 
low. Initial configurations were prepared by equilibrat- 
ing at least 1000 independent systems in the supercooled 
temperature regime, followed by quenching to the target 
temperature at a rate of lO -4 -^:. If not stated oth- 
erwise all the simulation below were obtained with sys- 
tems of N — 20164 in 2-dimensions and N = 16384 in 
3-dimensions. 

We choose to measure the sum 
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using an athermal, quasi-static scheme [10]. This scheme 
consists of imposing the affine transformation (4) to each 
particle of a minimized configuration, followed by an- 
other potential energy minimization under Lees-Edwards 
boundary condition [12]. In athermal quasi-static con- 
ditions (T — > 0, 7 — > 0), the system lives in local 
minima, and follows strain-induced changes of the po- 
tential energy surface [14]. Therefore, the particles do 
not follow homogeneously the macroscopic strain, and 
their positions change as — > r[ + Ui, where Ui de- 
notes non-affine displacements. Around some stable ref- 
erence state at 7 = 7*, the field it^, the potential energy 
U , and internal stress a xy are smooth functions of 7. 
We choose the stopping criterion for the minimizations 
to be |VjJ7| < 10~ 9 -| for every coordinate xt. Within 
this method one can obtain purely elastic trajectories of 
stress vs strain [10]. To measure B 2 of a given configu- 
ration of our molecular dynamics simulation at any tem- 
perature T, we first cool that configuration to T = 10~ 3 
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FIG. 3: Color Online: Upper panel: Trajectories of stress vs. 
strain for four different temperature at the same strain rate 
7 = 10 -4 . Lower panel: the corresponding values of B2 as 
a function of strain. Data was averaged over 1000 indepen- 
dent stress-strain curves at each temperature. Note that B2 
is negative even when the averaged stress-strain curve has a 
positive curvature, see text for discussion. 



using molecular dynamics during a time interval of 50r. 
This chosen temperature is sufficiently low to exclude 
any thermal activation on the time scale of the simula- 
tions. This initial treatment brings our configuration to 
an elastically stable state i.e. a minimum of the potential 
energy landscape. Without doing so one can find oneself 
in the vicinity of a saddle point for which the athermal 
elastic moduli have no clear meaning. We then apply the 
athermal quasi-static scheme to measure the finite differ- 
ences approximation to d £V ~ <Tx «^ 7 ) +o ' 3! ^y 7 ' > ~ 2 ' Tx!, ( ) 
by sampling a small elastic trajectory, using strain incre- 
ments of^7 = 2.5xl0~ 6 . We have checked that stricter 
stopping criteria for the minimizations or smaller strain 
increments do not significantly alter our results. We em- 
phasize that although we measure B2 in the athermal 
limit, the configurations on which we perform this mea- 
surement are sampled from various finite temperatures, 
see below. This athermal measurement is motivated by 
the requirement to probe the purely mechanical response, 
excluding thermal activation effects on the measurement 
from the discussion. Using this method we can compute 
B2 at any point of the trajectory. Note that B 2 is still a 
strong function of the temperature from which the con- 
figuration was taken, and this is because the organization 
of the particles depends on the temperature. We reiter- 
ate that B2 is not the second derivative of the averaged 
stress-strain curve, but rather the mechanical response of 
the underlying inherent structure which is sampled at a 
given temperature. 
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FIG. 4: Color Online: The Bauschinger effect for the four 
temperatures shown in Fig. 3, increasing from top to bottom. 
The trajectories are displaced by fixed amount (Alcr^l = 
0.15) for clarity. Note the reduction of the effect with increas- 
ing temperature. Data was averaged over 500 independent 
stress-strain curves at each temperature. Inset: the shaded 
area of difference between the stress-strain curves with posi- 
tive and negative strain as a function of B^- The magnitude 
of the Bauschinger effect saturates for T — I 0. 



III. RESULTS AND DISCUSSION 

The average trajectories of both stress vs strain (upper 
panel) and B2 vs strain (lower panel) are shown for four 
different bath temperatures in Fig. 3; the system was 
strained until 7=1/2 and then strain was reversed until 
the mean stress dropped to zero. The strain value was 
then 70 from which the experiment in Fig. 1 right panel 
was started with positive and negative straining with re- 
spect to 7o- The resulting trajectories of stress vs. strain 
are shown in Fig. 4 for the 2D system at the same four 
values of the temperature as in Fig. 3. We observe that 
the value of B 2 at the point of zero stress 70 reduces 
when the temperature increases, and in accordance with 
that the magnitude of the Bauschinger effect goes down 
as seen in Fig. 4. 

We can draw the conclusion that the magnitude of B 2 
is correlated with the amplitude of the Bauschinger effect 
(measured as the area of difference between the positive 
and negative stress-strain curves, see inset in Fig. 4). But 
even more detailed information which is highly relevant 
to the elasto-plastic behavior can be gleaned from the 
probability distributions functions (pdf's) of B 2 - These 
pdf's have rich dynamics along the stress-strain curves, 
as can be seen in Fig. 5. When measured in the isotropic 
zero-stress systems that are freshly quenched the distri- 
bution is symmetric as expected, with zero mean. In 
the elasto-plastic steady state the distribution moved to 
have a negative mean, in accordance with the low panel 
of Fig. 3. In Sect. IV we show that this distribution 
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FIG. 5: Color Online: Probability distribution function of B2 
at various points on the stress-strain trajectory. Data was 
collected from 3000 independent stress-strain trajectories at 
T = 0.01. In red (continuous) line we draw the symmetric 
pdf of the freshly prepared samples with 7 = 0. In green 
(dashed) line we show the pdf in the steady state, where it 
gains a negative asymmetry. In blue (dashed-dotted) line 
we see the pdf at the Bauschinger point 7 = 70 where it 
gained a positive asymmetry. The dynamics of these pdf's 
and their means are correlated with the shapes of the stress- 
strain curves and are proposed to be a crucial ingredient in 
any theory of elasto-plasticity. Inset: the N dependence of 
the pdf at the Bauschinger point 7 = 70. Data was averaged 
over 1000 independent samples for each system size. 



must send a tail towards —00 to accommodate the sharp 
changes in first derivative (the shear modulus) due to the 
proximity of mechanical instabilities in the form of plas- 
tic drops [13, 14]. At the Bauschinger point 70 the mean 
stress is zero, but the pdf of B2 gains a positive asym- 
metry, sending a tail towards +00, signalling a proximity 
to a plastic event in the negative straining direction. In 
the inset of Fig. 5 we exhibit the size dependence of the 
pdf at the Bauschinger point, to show that the asymme- 
try and the general shape of the pdf is quite independent 
of the number of particles N, always having long tails, 
indicating that near the Bauschinger point 70 there are 
close-by lurking plastic instabilities that are heralded by 
the tail of our pdf. 



To confirm that the qualitative findings reported above 
remain unchanged in 3-dimcnsions we repeated similar 
simulation for the model described above. In Fig. 6 we 
present a representative averaged stress-strained curve in 
the upper panel and the corresponding trajectory of B 2 , 
both at T = 0.01. 
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FIG. 6: A representative averaged stress-strain curve (aver- 
aged over 600 independent trajectories) for the 3-dimensional 
model (upper panel) and the corresponding trajectory of B2 
in the lower panel for T = 0.01. 



Finally, we use the crucial assumption [14] that the 
eigenvalue Ap crosses zero with a finite slope in the X- 
coordinate system itself, where distances are measured 
along the unstable direction: 



ApK AX + 0{X 2 



(13) 



Together with Eq. (12) and asserting that Sj is not sin- 
gular (it is a combination of derivatives of the potential 
function [17]), implies that 



X(-y) oc V7p - 7 



(14) 



These results are now used to determine the singularity 
of the stress at jp. We start with the exact result for the 
shear modulus [15, 16] 
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IV. THE DIVERGENCE OF B 2 NEAR 
MECHANICAL INSTABILITIES IN THE FORM 
OF PLASTIC EVENTS 

To see that B2 must reach ±00 when the system goes 
through a plastic deformation recall that as long as the 
system remains in athermal mechanical equilibrium (i.e. 
along the athermal elastic branch) the force /j on every 
particle is zero before and after an infinitesimal deforma- 
tion; In other words [14] with U the potential energy 



dfi _ d_dU_ _ d dU 
d'y d'y dTi d'y dui 
d 2 U d 2 U du. 
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where summation is implied by repeated indices. This 
condition introduces the all-important Hessian matrix 
Hij and the 'non-affine force' Sj which can both be com- 
puted from the interparticle interactions. We rewrite this 
condition as 
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where the second equation results from expanding in the 
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eigenfunctions of H, Hijipj = \kw\ 1 an d the last 
estimate stems from our knowledge that in finite systems 
the plastic event is associated with a single eigenvalue 
going through zero when the systems slides over a saddle. 
We denote the critical eigenvalue as Ap. Eq. (11) can be 
integrated to provide the distance of the non-affine field 
itj from its value at jp, 1^(7) — itj(7p) = X{^)ipf\ 
where X(j) is a function of 7 only, satisfying 



dX( n ) 
d'y 



Ap 



(12) 



where /ip is the Born term /ip = V~ 1 d 2 U/dj 2 . Using 
Eqs. (13) and (14) we conclude that near jp we can write 
the shear modulus as a sum of a regular and a singular 
term, 



a/2 



+ C(V7P ~ 7) 



(16) 



Obviously, our second derivative B2 = will inherit 

the singularity from the first derivative, explaining the 
long tails of the distributions seen in Fig. 5. Close to 
mechanical instabilities B2 is expected to diverge like 
(7p — 7)~ 3 ^ 2 with a sign that depends on the direction 
of imposed strain. 

Even though we sample our configurations from sys- 
tem with temperature where the singular points are not 
reached due to thermal activations, the proximity of these 
mechanical instabilities in a specific straining direction is 
signalled by the large values of B2 that we measure. 



V. SUMMARY AND CONCLUSIONS 

We have proposed here a new measure of the 
deformation-history induced anisotropy in amorphous 
solids. This measure is not model dependent and is eas- 
ily accessible to simulations and experiments. It is not 
obvious at this point in time whether a theory of elasto- 
plasticity should take into account the full pdf of B2, or 
whether it would be sufficient to take the mean value of 
B2 into account. We propose however that this object 
and its pdf are tempting analogues of the object m of 
the STZ theory as discussed above, with the obvious ad- 
vantage that they can be easily measured. In fact, in a 
follow up paper [17] we will show that this object can 
be expressed as a sum over the particles in the system, 
and therefore the measurements of the pdf can be done 
naturally and rapidly, making them highly accessible for 
further research. We stress that the value of B2 which 
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has been denned as the limit T — > in Eq. (9) can be 
measured experimentally at sufficiently low temperatures 
where the Bauschingcr effect is expected to be saturated. 
It appears worthwhile to measure this quantity in such 
low-temperature experiments and to correlate the value 
with the amplitude of the Bauschinger effect. 
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